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ABSTRACT 

We present detailed analyses of oxygen K absorption in the interstellar 
medium (ISM) using four high-resolution Chandra spectra towards the X-ray 
low-mass binary XTE J1817-330. The 11-25 A broadband is described with a 
simple absorption model that takes into account the pilcup effect and results in 
an estimate of the hydrogen column density. The oxygen K-edge region (21-25 A) 
is fitted with the physical warmabs model, which is based on a photoionization 
model grid generated with the xstar code with the most up-to-date atomic 
database. This approach allows a benchmark of the atomic data which involves 
wavelength shifts of both the K lines and photoionization cross sections in order 
to fit the observed spectra accurately. As a result we obtain: a column density 
of A h = 1.38 ± 0.01 x 10 21 cm -2 ; ionization parameter of log£ = —2.70 ± 0.023; 
oxygen abundance of Aq = 0.689igmoi and ionization fractions of O I/O = 
0.911, O Il/O = 0.077, and O III / O = 0.012 that are in good agreement with 
previous studies. Since the oxygen abundance in warmabs is given relative to 
the solar standard of Grevesse & Sauval (1998), a rescaling with the revision by 
Asplund et al. (2009) yields A 0 = 0.952+q;q 2 3, a value close to solar that rein- 
forces the new standard. We identify several atomic absorption lines — K cc, K f3, 
and Ky in O I and O II; and Kg: in O III, O VI, and O VII — last two probably 
residing in the neighborhood of the source rather than in the ISM. This is the first 
firm detection of oxygen K resonances with principal quantum numbers n > 2 
associated to ISM cold absorption. 

Subject headings: atomic data — atomic processes — ISM: general — X-rays: general 
— X-rays: binaries — stars: individual (XTE J1817-330) 
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1. Introduction 

X-ray studies of interstellar absorption have the potential to provide information that 
is not accessible by other techniques. This includes the relative abundances of a wide range 
of ion stages of interstellar cosmic elements (e.g. oxygen) and the abundances relative 
to H and He of elements with atomic number Z > 6. X-ray absorption can also provide 
signatures of the binding of these elements in molecules or solids; the inner-shell electronic 
transitions are key diagnostics since the ionization state or chemical binding shifts the line 
energies by a predictable amount. The ongoing effort to explore and exploit this effect in 
interstellar studies involves accumulation of synthetic spectra for inner-shell absorption of 
ions and neutrals, based on both atomic calculation and laboratory experiment, in addition 
to testing them against available observed interstellar X-ray absorption spectra. 

As part of the effort to accumulate the required atomic data, the atomic database of the 
xstar photoionization modeling code (Bautista & Kallman 2001) has been systematically 
improved in the past ten years in order to study the K lines and edges in the high-quality 
X-ray astronomical spectra obtained from the Chandra and XMM-Newton satellite-born 
observatories. This database includes K-vacancy levels, wavelengths, A-values, radiative 
and Auger widths and high-energy photoionization cross sections for complete isonuclear 
sequences with Z < 30 (see Palmeri et al. 2012, and references therein). 

With this new database it has been possible, for instance, to estimate the model 
efficiency of iron K-line emission and absorption in terms of the ionization parameter and 
density (Kallman et al. 2004). It has been shown that the centroid of the Ka unresolved 
transition array, the K/3 energy, and the ratio of the Kaq to K «2 components are useful 
diagnostics of the ionization parameter. It was also found that the many strongly damped 
resonances below the K-ionization thresholds lead to edge smearing, which can certainly 
hamper the astrophysical interpretation of the absorption features. A synthetic spectrum 
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has been modeled by Kallman et al. (2009) for the black-hole X-ray transient GRO J1655-40 
which constrains the odd -Z element abundances (11 < Z < 27) and outflow parameters 
of the associated warm absorber. By analyzing the oxygen K-absorption structure in the 
interstellar medium (ISM) from the XMM-Newton spectrum of the low-mass X-ray binary 
Sco X-l, Garcia et al. (2011) reproduced both the K edge and O I Ka absorption line, 
thus evaluating the impact of the atomic data on the interpretation of observations. In the 
present report we continue this work with a more stringent benchmark of the XSTAR atomic 
database, made possible by Chandra spectra of a low-mass binary displaying imprints of 
the ISM oxygen K absorption with higher resolution. 

Juett et al. (2004) measured high-resolution spectra of the interstellar oxygen K-shcll 
absorption edge from seven X-ray binaries — Cygnus X-l, Cygnus X-2, 4U 1636-53, 

4U 1735-44, GX 9+9, 4U 1543-624, and 4U 1820-30 — using the Chandra High Energy 
Transmission Grating Spectrometer (HETGS). The oxygen column density Nq was therein 
calculated using the optical depth at 21.7 A and the cross section at this wavelength 
assuming the oxygen abundance from Wilms et al. (2000). With this value, the hydrogen 
column density iV n was then determined. The Is- 2p transitions of O II and O III were 
identified, finding ISM abundances for O II and O III relative to O I of ss 0.1 and <0.1, 
respectively. This became the first estimate of the O II/O I and O III/O I abundance ratios 
in the ISM. 

By observing the oxygen K-shell edge toward galaxy clusters, Baumgartner & Mushotzky 
(2006) showed that the ISM hydrogen X-ray column density was in close agreement with 
the 21 cm radio value for columns less than approximately 0.5 x 10 21 cm -2 . For higher 
column densities, the former was higher by as much as a factor of 2.5, indicating substantial 
absorption besides that due to neutral hydrogen (probably from clouds of molecular 
hydrogen). An average ISM oxygen abundance of 0.99 solar was found relative to the 
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photospheric value (8.66 ± 0.05) by Asplund et al. (2004) which implied a high gas to dust 
ratio. 

Pinto et al. (2010) reported an spectral model suggesting that, along the line of sight 
towards an X-ray source (GS 1826-238), the ISM was more complex than just simple neutral 
gas. Their gas model consisted of three components: cold gas (T « 5.8—10 x 10 3 K); warm 
ionized gas (T « 1—6 x 10 4 K); and hot ionized gas (T « 2 x 10 6 K). It was also shown 
that the column densities N h of these three components span two orders of magnitude, 
where the cold gas contributed ~ 90—95%, the warm gas ~ 5—10%, and the hot gas ~ 1% 
to the total N^. Costantini et al. (2012) analyzed the ISM using the bright X-ray binary 
4U 1820-30 source, concentrating on the mildly ionized cold gas where the model consisted 
of mildly ionized gas, dust, and mildly ionized outflowing gas. In the oxygen edge region, 
dust compounds were included which did not show features below 23.7 A. Finally, oxygen 
was found to be overabundant by a factor of 1.23 solar. 

In order to model the oxygen K absorption features in the ISM with xstar, we need a 
Chandra spectrum with high counts. In the catalog of low-mass X-ray binaries by Liu et al. 
(2007), we found in XTE J181 7-330 a suitable source with an average count number of 
~ 8 x 10 6 counts/s. Discovered by the Rossi X-ray Timing Experiment (Remillard et al. 
2006), XTE J1817-330 is a low-mass X-ray binary displaying a primary component which 
has an X-ray spectrum and variability behavior which are consistent with its identification 
as a black hole (Sala et al. 2007). Its location, l = 359.81 and b = —7.99 in Galactic 
coordinates, implies a line of sight in the Galactic center direction, a region that dominates 
the total mass of the Galaxy (Stasinska et al. 2012) providing a good opportunity to 
study the ISM absorption features. This source has been observed in different spectral 
ranges, namely the optical (Torres et al. 2006), near infrared (D’Avanzo et al. 2006), 
radio (Rupen et al. 2006), and X-ray (Roy et al. 2011). Sala et al. (2007) fitted its X-ray 
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spectrum finding an absorbing hydrogen column of 1.55 x 10 21 cm -2 and identifying the 
following oxygen absorption lines: O I Kct(23.52 A); O II Kct(23.35 A); O III Ka/23.13 A); 
O I K/3(22.91 A); and O VII K/3(21.60 A). Gierlinski et al. (2009) estimated a hydrogen 
column density of X H = 1-08 ± 0.05 x 10 21 cm" 2 , and Roy et al. (2011) investigated the 
timing characteristics in the spectra using a column of Ah = 1.2 x 10 21 cm" 2 . 

The outline of the present report is as follows. In Section 2, we describe the xstar 
warraabs physical model which is used to obtain the final spectrum fits, in particular the 
details concerning its atomic database. The data reduction aspects of the observations, 
e.g. the pileup effect, are reviewed in Section 3 while the results of the broadband fits and 
resolved absorption features are given in Sections 4-5. Section 6 is dedicated to comparisons 
of the hydrogen column density that allow an estimate of its reliability. Section 7 describes 
the oxygen column densities derived from our fits. Finally, a discussion and conclusions are 
presented in Section 8. 


2. XSTAR/ warmabs 

xstar 1 is a code package designed to determine the physical conditions in partially 
ionized gases. It calculates the temperature, ionization, excitation, absorption, and emission 
due to neutral and all ionized species of elements with atomic number Z < 30. As many 
relevant physical processes as possible are taken into account, assuming a time stationary 
balance among them and a Maxwellian electron velocity distribution. Important processes 
for many problems relevant to X-ray astronomy include photoionization, electron impact 
collisional ionization and excitation, and radiative and dielectronic recombination. More 
complete descriptions of xstar and its atomic database are given by Kallman & Bautista 


1 http://heasarc. nasa.gov/lheasoft/xstar/xstar.html 


(2001) and Bautista & Kallman (2001). In this work we employ the warmabs utility 
which allows for the calculation of XSTAR models within the widely used X-ray spectral 
fitting package xspec 2 . It utilizes precalculated tabulations of the ion fractions and level 
populations in order to calculate synthetic emission and transmission spectra. It assumes 
that the gas responsible for emission or absorption has a uniform ionization and temperature 
throughout, although it is possible to superimpose two or more components and simulate a 
non-uniform situation. 


2.1. Atomic data 

The K photoabsorption data sets for oxygen ions with electron numbers 2 < N < 8 
contained in the xstar atomic database were generated by Garcia et al. (2005). Energies 
for both valence and K-vacancy fine-structure levels, transition wavelengths, (//-values, 
and radiative and Auger widths were computed with the atomic structure packages 
AUTOSTRUCTURE (Eissner et al. 1974; Badncll 1986, 1997) and hfr (Cowan 1981). 
Electron correlation effects were taken into account by means of configuration-interaction 
expansions within the n < 3 complexes, and relativistic corrections were included with a 
Breit-Pauli Hamiltonian. The accuracy of the theoretical level energies and wavelengths 
does not generally match that attained in experiments (see Table 1), and the scanty 
availability of measurements for the K shell of O ions with N > 3 limits the implementation 
of term-energy corrections and wavelength shifts that can be introduced in order to 
compensate for the theoretical shortcomings. Furthermore, as also shown in Table 1, 
experimental misassignments, such as that quoted by Gu et al. (2005) for O VI, are not 
uncommon. (The discrepancy in O III is believed to be a typo.) 


2 http: / / heasarc.gsfc.nasa.gov/xanadu/xspec/ 
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High-energy photoionization cross sections in the K-edge region for oxygen ions 
were obtained by Garcia et al. (2005) with bprm, the Breit-Pauli A-matrix suite of 
codes (Scott & Burke 1980; Scott & Taylor 1982; Berrington et al. 1987; Seaton 1987). A 
conspicuous effect in the resonance structure of these cross sections that must be treated 
in detail is Auger damping. It is the result of the dominant spectator KLL Auger decay 
channels of a photoexcited K-vacancy state, say, 


[ls]2pGip — » 

2p M_2 np + e" 

(1) 

->■ 

[2s]2p M-1 np + e - 
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which leads to edge smearing by resonances with symmetric profiles of nearly constant 
width. They are handled within the A-matrix package by means of an optical potential 
(Gorczyca & Robicheaux 1999; Gorczyca & McLaughlin 2000). Again, the accuracy of the 
resonance energy positions is limited by the absence of measurements which would enable 
empirical adjustments of the many K-vacancy target thresholds (i.e. the series limits) in 
the close-coupling expansion. 


3. Observations and data reduction 

In Table 2 we list the specifications of the four observations of XTE J1817-330 obtained 
by Chandra using the High Energy Transmission Grating Spectrometer (HETGS) in 
combination with the Advanced CCD Imaging Spectrometer (ACIS). The oxygen edge is 
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accessible with the HETGS exclusively via the Medium Energy Gratings (MEG). Given 
the large flux of this source (F = 1350 /i J, integrated in the 2-10 keV energy range), 
these observations were taken in continuous clocking mode which significantly increases the 
temporal resolution in order to minimize the pileup effect (Cackett et al. 2008). The pileup 
is an inherent feature of CCD detectors, such as those of the ACIS instrument, which occurs 
when two or more photons are detected as a single event possibly causing a deformation in 
the level and shape of the continuum (Miller et al. 2006b). It is stronger in the MEG than 
in the HEG due to the lower dispersion and higher effective area of the former 3 . Although 
the pileup is usually not an issue near the oxygen edge, it may be for the continuum which 
is used to establish the hydrogen column density by fitting the 11-25 A interval; hence, we 
consider the pileup effect in this range by applying the correction model simple_gpile2 . si 
(Hanke et al. 2009) before spectral fitting. Among the most popular computer packages 
for fitting X-ray data are xspec and ISIS, and we have used the latter (version 1.6.2 4 ) to 
include the pileup model. 

Even though the spatial resolution is reduced to one dimension when the continuous 
clocking mode is used, there is essentially no difference with the timed exposure mode in 
the background extraction or in the data reduction process. Hence, we have reduced the 
data sets using the standard CIAO threads 5 . In some cases, the zero-order data were not 
telemetered; for these observations the zero-order position was then estimated by finding 
the intersection of the grating arms using the f indzo algorithm 6 . 

Figure 1 shows the light curves of the four observations ObsID 6615, 6616, 6617, and 

3 http: / / cxc.harvard.edu/ciao/ ahelp/acis_pileup.html 

4 http:/ /space. mit.edu/cxc/isis/ 

5 http: //cxc. harvard. edu/ciao4.4/threads/gspec.html 

6 http : / / space . mit . edu / cxc / analysis / findzo / 
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6618 in units of counts/s. It may be seen that the average counts/s is approximately 
constant for all, and therefore, we use the averaged spectra for this work. We have carried 
out the following spectral fitting procedure: firstly, each spectrum was rebinned with a 
0.1 A bin in order to fit the continuum in the 11-25 A wavelength range so as to find the 
hydrogen column density using the TBnew model; then, the default spectral resolution was 
restored, and the oxygen edge region (21-24 A) was analyzed with the warraabs model 
freezing the column density at the value obtained in the preceding step. 


4. Broadband fit 

Following Yao & Wang (2006), we rebinned each spectrum with a 0.1 A bin size 
throughout the spectral range to decrease the number of absorption lines and improve 
the broadband fit. We fitted the four observations simultaneously in the interval 11-25 A 
using the simplegpile2 (TBnew (powerlaw) ) model, where simplegpile2 is a convolution 
model to account for pilcup effects and TBnew is an X-ray absorption model that includes 
abundances for elements from H to Ni (Wilms et al. 2000). Also we added Gaussians to fit 
the remaining absorption lines after rebinning. We used abundances specified in Wilms et al. 
(2000) and cross sections by Verner et al. (1996). The O, Ne, and Fe abundances were 
handled as free parameters, and for the analysis \ 2 statistics was employed. We fixed 
the absorption parameters in all the observations, and varied the power-law, pileup, and 
Gaussian parameters for each case. 

Figure 2 shows the fit in the interval 11-25 A using the simplegpile2(TBnew(powerlaw) ) 
model. Although the fit is carried out simultaneously for the four observations, they 
are plotted separately to improve clarity. Panels (a), (b), (c), and (d) correspond to 
observations ObsID 6615, 6616, 6617, and 6618, respectively (see Table 2). In each panel, 
black data points correspond to the observation while the solid red lines represent the 
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best-fit model for each case, and the base plots show the fit residuals in units of y 2 . In each 
spectrum the K edges of Ne and O are respectively located at ~ 14.3 A and ~ 23.2 A. The 
O I and O 11 Is- 2p absorption lines at 23.5 A and 23.35 A are also clearly observed. We 
improved the fit by including Gaussians to model the outstanding absorption lines: Ne IX 
Kcc (« 13.43 A); Ne III Kct (« 14.50 A); Ne II Ka (« 14.60 A); O VIII Lycc (« 18.95 A); 
O VII K a (~ 21.58 A); and O II Ka (~ 23.34 A). Residuals indicate that the model fits 
the data satisfactorily with deviations near the Ne edge (~ 14.3 A) due to background (i.e. 
clocking the whole chip) or an unaccounted hot pixel. We also find high residuals in the 
oxygen absorption region; in particular, the O I K a and O II Kct are not well modeled. 

Fit results for the individual free abundances, column density, power law, and pileup 
are presented in Table 3. The fit statistics give a reduced chi-square of y 2 = 1.208. 
Differences in the power-law parameters are due to the count numbers in each observation. 
We obtain abundances for O, Ne, and Fe higher than solar (Wilms et al. 2000), especially 
for Ne (factor of 2), and a column density of Nn = I.66I004 x 10 21 cm 2 which is somewhat 
larger than those estimated in other spectral fits and by measurements of the 21 cm line. 
This will be further discussed in Section 6. 

The nonlinear pileup convolution model siraple_gpile2 . si exponentially reduces the 
predicted count rate Rj( A) (in units of counts/s/A) according to 

R(\) = R(A)exp(-/TR tot (A)) (6) 

where R to t (A) is the total spectral count rate and f3 is treated as a fit parameter. The model 
is only applicable to first order grating spectra although it employs information from the 
second and third order grating spectra 7 . With this model we obtained an average (3 = 0.050 
according to which the spectra have a pileup degree greater than 25% for 14 < A < 16 A 


' http: / / cxc.harvard.edu/ciao/ahelp/acis_pilcup.html 
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with the highest value of p — 42% at ~ 14.3 A for the observation with the largest number 
of counts (ObsID 6615). This value shows the need to include the pileup effect in the 
11-25 A region to obtain a good fit. However, for all the observations, we have a pileup 
degree lower than 5% in the 21-24 A wavelength range which corresponds to the oxygen 
absorption region and, therefore, its effects were ignored. 


5. Absorption features 

In order to identify all the absorption features present in the 21-25 A wavelength 
region, we first applied a functional model consisting of a power law for the continuum and 
Gaussian profiles to describe the absorption lines. Figure 3 shows the spectral fit of the four 
Chandra MEG observations of XTE J1817-330 in the oxygen absorption region. Although 
shown in separate panels for clarity, the fit was performed simultaneously for the four 
observations. Nevertheless, the power-law parameters (photon index and normalization) 
are allowed to vary independently for each data set. Several absorption features are clearly 
observed, the most outstanding being the O 1 and O II Ka lines. The base plot in each 
panel shows in units of y 2 the fit residuals as a histogram. High residuals are found around 
24.8 A, and the model underestimates the number of photons around 23.2 A. 

Table 4 shows a list of the lines obtained with the functional fit including the centroid 
wavelengths derived from the Gaussian profile parameters. It also lists the wavelengths 
for the lines measured by Sala et al. (2007). Note that our wavelengths are within the 
error bars of the values reported by these authors, which are much larger given the lower 
spectral resolution of the XMM-Newton instruments. We obtain a reduced chi-square 
of y 2 = 1-434. The O I Ka is found at 23.502 ± 0.001 A which has been previously 
declared of ISM origin (Sala et al. 2007). Furthermore, we have also detected four other 
O I absorption lines corresponding to resonances with higher principal quantum number n: 
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two pairs of K/3 (n = 3) and K7 (n = 4) resonances associated to the 4 P and 2 P K-hole 
core states. These are found at 22.884 ± 0.004 A, 22.790 ± 0.001 A, 22.686 ± 0.004 A, and 
22.609 ± 0.007 A, respectively. In this set, only the first K/3 (belonging to the 4 P core state) 
has been previously detected at 22.91 ± 0.03 A in the XMM-Newton observation analyzed 
by Sala et ah (2007). This is then the first firm detection of oxygen K resonances with 
n > 2 associated to ISM cold absorption. 

The line at 23.358 ± 0.002 A is also of ISM origin (Juett et ah 2004) and corresponds 
to O II K a. As in the case of O I, resonances with n > 2 were also detected for O II. 
These were found at 22.280 ± 0.003 A and 22.101 ± 0.005 A corresponding to K/3 and Ky, 
respectively. Although the O III K a resonance is a triplet, only two absorption features 
were detected at 23.104 ± 0.005 A and 23.054 ± 0.001 A while Sala et ah (2007) only 
found one line at 23.13 ± 0.09 A identified as O III Ka. Finally, we have also located Ka 
transitions for both O VI (22.022 ± 0.003 A) and O VII (21.589 ± 0.003 A), from which only 
the strongest of the two (O VII) was reported by Sala et al. (2007) at 21.609 ± 0.06 A. 

The functional fit described above illustrates the complexity of oxygen K absorption. 
There are (at least) two different regimes: a cold component of mostly neutral gas and a hot 
component with highly ionized gas. Previous studies have suggested a multiple-phase ISM; 
however, it is not clear if the two high-ionization lines from O VI and O VII are intrinsic 
to the source or in fact arise in the ISM. Thus we direct our study to the cold phase. 
Furthermore, since cold absorption is not exclusively due to neutral oxygen, this means 
that a model such as TBnew will not completely represent the observed features. We have 
then simultaneously fitted the four exposures in the 21-25 A wavelength range using the 
powerlaw*warmabs physical model. The advantage of using warmabs is that it contains the 
most recent atomic data for all the ions in the oxygen isonuclear sequence, thus enabling 
the fit of all the absorption lines from O I, O II, and O III. While implementing this model, 
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the column density is held fixed at the value obtained in the broadband simultaneous fit 
(A r H = 1.66 x 10 21 cm -2 ). The ionization parameter and the oxygen abundance are treated 
as free parameters but the same for all observations, while the abundances of all the other 
chemical elements are held fixed at their solar values. The quantities describing the power 
law (photon index and normalization) are also taken as free parameters but independently 
for each observation. 

Figure 4 shows the results of this fit. As before, we present for clarity the spectra 
with the best-fit model in separate panels although the fit is performed simultaneously. 
Since most of the relevant absorption features in the spectra are due to neutral, singly, and 
doubly ionized O ions, the ionization parameter that describes the fit tends to be somewhat 
low: log £ = —2.729. In this case, the warmabs model does not show any absorption due to 
O VII; therefore, this line is still represented with a Gaussian profile. We have found large 
residuals around most of the oxygen absorption features, the most prominent being around 
the O I and O II Kcc lines. The fit statistics gives a poor reduced chi-square of y 2 = 1.769 
which forced us to revise the atomic database. 

A close examination of the resonance positions from the photoabsorption cross sections 
used in warmabs in comparison with those obtained from the spectral fits with Gaussian 
profiles reveals the inaccuracies of the atomic data. Figure 5a shows with solid colored 
lines the atomic cross sections used in warmabs for O I, O II, and O III (Garcia et al. 
2005). Vertical dashed lines are placed at the positions of the absorption features found 
with the functional fit described above (see Table 4). The comparison shows that not 
only the positions of the O I and O II K a resonances are displaced with respect to the 
observed data but also those with n > 2. The solid black line represents the experimental 
O I photoabsorption cross section measured by Stolte et al. (1997) (unfortunately, there 
are no laboratory measurements available for the other species) which displays a similar 
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wavelength shift with respect to the observed resonance positions. This wavelength offset 
has been previously discussed by Juett et ah (2004) based on the Ka position observed in 
the Chandra spectra from several sources, all of which match the position we have found 
using Gaussian profiles. For O I, the comparison of the theoretical cross section with the 
laboratory measurement shows that, except the Ka peak, all the other resonances agree 
very well. This implies that the relative position of the Ka with respect to the other 
resonances also needs to be adjusted. This is also the case for O II since the wavelength 
difference of the theoretical and observed Ka positions is not the same as those for the 
rest of the resonances of the same ion. Consequently, we have decided on the one hand to 
adjust the positions of the Ka resonances (which are treated as lines in warmabs), and on 
the other, to shift the whole cross sections for both species so as to obtain the best possible 
agreement with the observed lines. The new positions for the O I and O II Ka lines are 
23.502 A and 23.343 A while the cross sections were shifted towards shorter wavelengths 
by 0.033 A and 0.079 A, respectively. No correction was applied to the O III curve since 
the observed lines are weaker, in addition to the fact that in this case the Ka resonance 
is a triplet. Figure 5b shows the atomic absorption cross sections after the wavelength 
corrections are applied. In other words, the Chandra observed line positions suggest not 
only that the theoretical resonance positions from Garcia et al. (2005) must be adjusted, 
but also that the wavelength scale of the only available experiment for O I (Stolte et al. 
1997) needs to be shifted by 33 mA. 

Possible effects due to the instrument wavelength calibration, model uncertainties, 
and Doppler shifts due to motion of the gas in the ISM are important for accurate 
line identification. The instrumental resolution for the first-order spectra from MEG is 
AA ~ 23 mA, which is better than the largest shift required for our theoretical cross 
sections (33 mA). In fact, according to Garcia et al. (2005), the uncertainties in the 
theoretical cross sections for O ions with electron occupancies N < 4 were estimated to be 
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~ 50 mA, comparable with the wavelength uncertainties in the laboratory measurements 
of Stolte et al. (1997). Assuming a velocity dispersion for the ISM of < 200 km s -1 
(Juett et al. 2004), the largest wavelength shift due to Doppler effects for the O I Kct is 
AA = 14 mA, which is smaller than the instrumental resolution. Hence, we can safely 
ignore the effects of ISM gas motion in the light of sight. 

Figure 6 shows the fit results with warmabs after the theoretical line positions have been 
corrected, where the residuals are now significantly reduced and more evenly distributed. 
This is reflected in better fit statistics (reduced y 2 = 1.245) and a well-constrained ionization 
parameter (log£ = —2.699 ± 0.023) which is consistent with the low-ionization oxygen lines 
observed. As before, the O VII Ka resonance at 21.589 A is not well represented with this 
model and must be fitted with a Gaussian profile. We have also noticed that the absorption 
feature at 22.022 A is not entirely due to the O III K<5 resonance and is thus fitted with a 
Gaussian profile, confirming a blend of both O III KA and O VI Kq. In this final analysis, 
we also treated the column density value as free parameter in order to increase the accuracy 
of the fit. We find a column density of Ah = 1.38 ± 0.1 x 10 21 cm -2 that is lower than 
the value previously mentioned in connection with the broadband fit (see Section 4). The 
oxygen abimdance obtained with the warmabs model is Aq = 0.689^q q(q relative to solar, 
also lower than the value for the best fit with TBnew (Ao = 1.178 ± 0.22). All the best-fit 
parameters of the warmabs are summarized in Table 5. 


6. Hydrogen column density 

Table 6 shows a comparison of the hydrogen column density obtained in this work 
with previous estimates including measurements using the 21 cm line. By means of the 
simultaneous broadband fit described in Section 4, we have obtained a column density of 
A h = 1.661q;q4 x 10 21 cm -2 which is somewhat larger than those quoted in previous spectral 
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fits. Miller et al. (2006a) reported a hydrogen column density of IVh = 8.8— 9.7 x 10 20 cm -2 
using a 50 ksec Chandra spectrum, a value smaller than those obtained from Galactic 
H 1 surveys, e.g. 1.58 x 10 21 cm -2 (Dickey & Lockman 1990) and 1.39 x 10 21 cm -2 
(Kalberla et al. 2005). Using the equivalent widths of several interstellar bands in the 
optical spectrum, Torres et al. (2006) found the column density to be in the range 
IVh = 1—3 x 10 21 cm -2 . Moreover, both Rykoff et al. (2007) and Gierlinski et al. (2008) 
gathered acceptable fits of several Swift XRT spectra with a fixed Galactic column density 
of IVh = 1.2 x 10 21 cm -2 . Gierlinski et al. (2009) used the same XRT spectra combined 
with data from UVOT (The Ultraviolet and Optical Telescope) to derive a column density 
of IVh = 1-08 ± 0.05 x 10 21 cm -2 from broadband spectral fits. 

The agreement between TBnew and Sala et al. (2007) is around 7% despite the fact that 
their fits include data from XMM-Newton EPIC-Pn (0.6—10.0 keV), RGS1 (0.3— 2.0 keV), 
and OM covering a much wider energy range; furthermore, they model the foreground 
absorption with TBabs, an older version of TBnew (Wilms et al. 2000). On the other hand, 
the present fits of the oxygen K region (21—25 A) with the warmabs physical model yield a 
lower column density (IVh = 1-38 ±0.1 x 10 21 cm" 2 ) than that derived from our broadband 
fit and TBnew. This discrepancy is expected due to the different atomic data sets involved, 
particularly since TBnew only includes the photoabsorption cross section for neutral oxygen. 
Taking into account the scatter of the previous column densities in the literature and the 
fact that we have adequately modeled the absorption features from O I, O II, and O III, we 
would expect the value derived from warmabs to be more reliable than that from TBnew. 


7. Oxygen column density 

In order to derive the O column density along the line of sight of XTE J1817-330, 
we have calculated the curve of growth for each of the oxygen Ka transitions using the 
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atomic data from Garcia et al. (2005), and equivalent widths (EWs) have been obtained 
using the functional fit described in Section 5. Table 7 shows a comparison between the 
present Ka EWs for O I, O II, O III, O VI, and O VII and those by Sala et ah (2007) for 
XTE J1817-330. EWs listed in Juett et ah (2004) for seven other sources are also included 
as well as those by Yao et ah (2009) for Cygnus X-2. For O I Ka, our EW (51 ± 5 mA) 
is around the lower limits of Sala et ah (2007) and Juett et ah (2004), and for O II Ka, 
the present value of 52 ± 5 mA is within the range that Juett et ah (2004) associated with 
the ISM. However, we End the Ka lines in both O I and O II to be saturated which could 
influence our derived EWs. Furthermore, the Ka EWs in O III, O VI, and O VII are in good 
agreement with the values obtained by Juett et ah (2004), Sala et ah (2007), and Yao et ah 
(2009). 

Table 8 shows a comparison of the oxygen column densities obtained from the EWs and 
the warmabs model fit. As previously mentioned in Section 5, for the ionization parameter 
found in our best warmabs fit, only O I, O II, and O III are included in the model. As 
a reference, we have also included the oxygen column density for O VI from Savage et ah 
(2003), obtained from FUSE data for four sources located near the same Galactic latitude as 
XTE J1817-330. The curve-of-growth values for O I, O II, and O III result in a total column 
density of No = Noi + Non + Nom = 5.71 x 10 1 ' cm -2 , in relatively good agreement with 
that derived from the ionic fractions in the warmabs fit, namely Nq = 6.41 x 10 1 ' cm -2 . 
However, the ratios O I/O II are significantly different: 1.4 and 11.9, respectively. Due to 
line saturation, the measured EWs are likely to be imderestimated, in particular for the 
O I Ka line. This leads to a higher uncertainty in the column and ion fractions estimated 
from the EW measurements. On the other hand, since the parameters in the warmabs fit 
are constrained by all the lines as well as the K edge, we are more confident about their 
reliability. In the case of the O VI and O VII Ka transitions, because the EWs are on 
the flat section of the curve of growth, they strongly depends on velocity dispersion. For 
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O VI Kcq the oxygen column density spans from Nq v1 = 71.12 ± 29.15 x 10 15 cm -2 to 
iVoyi = 5.06 ± 2.53 x 10 15 cm -2 when using dispersion velocities of v = 20 km s" 1 and 
v = 200 km s _1 , respectively. The smallest value is still about one order of magnitude 
larger than that reported by Savage et al. (2003), making difficult to constrain the the 
velocity dispersion. Also, no real improvement is achieved by considering larger dispersions, 
as the column density is reduced very little once the velocities are larger than 200 km s _1 . 
This may be an indication that the observed O VI line is likely to originate from the 
neighborhood of XTE J1817-330 rather than from the ISM. In the case of O VII, the 
oxygen column density varies from 7Vo vn = 4.48 ± 1.84 x 10 18 cm -2 (w = 20 km s^ 1 ) to 
N Qvn = 4.5 ± 1.8 x 10 16 cm- 2 (u = 200 km s" 1 ). 

8. Discussion and conclusions 

There are several interesting conclusions emerging from this work that configure a set 
of useful guidelines in the analysis of ISM absorption features in Chandra and XMM-Newton 
spectra. In order to resolve such spectral signatures, bright sources with high count rates 
are desirable. In the present XTE J1817-330 study case, we were fortunate to find four 
good quality spectra that were fitted simultaneously so as to obtain improved statistics. 
However, it was realized that the pileup effect must be taken into account in order to obtain 
reasonably reliable hydrogen column densities. This effect was treated adequately with the 
nonlinear pileup convolution model simple_gpile2 . si, whereby the pileup was estimated 
by fitting the continuum to be on average around 25%. The resulting hydrogen column 
density is V H = lffib^or x 10 21 cm 2 , which is somewhat (< 7%) higher than those by 
Kalberla et al. (2005) and Sala et al. (2007), the former obtained from 21 cm data. 

The absorption features in the oxygen region (21-25 A) were first picked up with a 
functional model (power law and Gaussian profiles) which yielded Ka lines of O I, O II, 
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O III, O VI, and O VII previously observed by Sala et al. (2007), but for the first time, 
also K resonances of O I and O II with principal quantum number n > 2. These line 
identifications at least confirm a two-phase plasma with both neutral (mostly) and highly 
ionized components. In the present work we have been mainly concerned with the former. 
In this respect, since neutral oxygen is not the sole charge state responsible for the cold 
oxygen absorption, the four exposures are then simultaneously fitted with a representative 
physical model referred to as warmabs, which contains the atomic data for the oxygen 
isonuclear sequence computed by Garcia et ah (2005). In this step, the hydrogen column 
density is held fixed at the previously determined value of A H = 1.66 x 10 21 cm -2 while 
the ionization parameters and oxygen abundance are treated as free parameters. A detailed 
comparison of this model with observations brings forth the inaccuracies of the atomic 
data, which led us to introduce wavelengths shifts in both the K a line positions and 
photoabsorption cross sections of O I and O II. 

The scope of this corrective procedure is severely limited by the general poor 
availability of laboratory data for the oxygen inner-shell, but it certainly leads to 
improved fits and statistics; in particular, to a well-constrained ionization parameter of 
log£ = —2.699 ± 0.023 consistent with the absorption lines observed. In order to further 
improve the fit, the column density is allowed to vary resulting in a more reliable value 
(Ah = 1-38 ±0.1 x 10 21 cm -2 ) which is within 13% of previous estimates (Kalberla et ah 
2005; Sala et al. 2007). Furthermore, the fit yields an oxygen abundance of A 0 = 0.689^q;q^q 
and ionization fractions of O I/O = 0.911, O II/O = 0.077 and O III / O = 0.012. Regarding 
the oxygen abundance, it may be clarified that the quoted value is given relative to 
the solar standard (8.83 ± 0.06) of Grevesse & Sauval (1998). If it is rescaled to the 
recently revised standard by Asplund et al. (2009) of 8.69 ± 0.05, then an ISM abundance 
(Ao = 0.9521S) close to solar is obtained that is in excellent accord with the averaged 
value by Baumgartner & Mushotzky (2006), and can be included as further support for the 
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solar abundance revision. Moreover, our ionization fractions are in good agreement with 
those found by Juett et ah (2004) of O II/O I ~ 0.1 and O III/ 0 I < 0.1. 

The functional fit allows us to derive oxygen column densities by means of a curve of 
growth method using the EWs for each of the observed Kct absorption lines, which can 
then be compared with the values obtained from the warraabs fit. Even though we find by 
both methods a relatively consistent total column density (5.71 x 10 17 cm -2 in the curve of 
growth and N = 6.41 x 10 17 cm -2 in the warmabs model fit), the O I/O II ratio for the Ka 
lines is significantly lower (an order of magnitude) in the former. This is attributed to the 
high saturation level in both lines which leads to underestimates of the measured EWs, and 
we are therefore confident that the results from the warmabs model fit are more reliable. 

In the case of O VI, the width of the line strongly depends on the velocity dispersion of 
the gas. For the velocity dispersion values expected in the ISM {y ~ 20 — 200 km s -1 ), 
the EWs obtained with our functional fit yield oxygen column densities far above from the 
FUSE values by Savage et al. (2003). This may be an indication that the O VI observed 
absorption in the XTE J1817-330 spectra is bound to occur near the source rather than in 
the ISM. 

We find acceptable fits to the observed oxygen K-edge spectrum with a model consisting 
of atomic oxygen ions if the energies are slightly adjusted. However, it is also of great 
interest to detect or set limits on features due to other phases of interstellar oxygen and 
other elements, namely molecules or solids (dust). It is expected that these phases comprise 
a significant fraction of the total oxygen column in the ISM. Evidence for this includes 
observed depletions of atomic oxygen along many lines of sight relative to certain elements 
(Jenkins 2009; Yao et al. 2009; Whittet 2010), together with observations of molecular 
features in other wavelength bands. An example of a potentially detectable signature of such 
material is the K resonance of the CO molecule expected at a rest wavelength of 23.21 A 
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(Barrus et al. 1979). This line is separated from the nearest feature we clearly detect by 
approximately 100 mA, which can be compared with the magnitude of the wavelength 
shifts we adopt in order to get acceptable fits, the largest of which is 33 mA. Therefore, we 
conclude that the existing atomic absorption data are not sufficiently uncertain to mask the 
presence of this molecular signature. The search for molecular and solid features in X-ray 
spectra remains of great interest, both as tests of the atomic data and of the models for the 
ISM. 

We believe that the present benchmark has produced an improved version of the 
warraabs model which can now be used for a more extensive and reliable study of oxygen 
K photoabsorption towards other X-ray sources. This will be the next step of the current 
study where the oxygen abundance variations in the ISM are certainly of interest. In a 
similar fashion, we also intend to evaluate the Ne and Fe edge regions for ongoing physical 
model refinement. 

Part of this work was carried out by Efrafn Gatuzz during attendance in July 2011 
to the Committee On Space Research (COSPAR) Capacity Building Workshop in San 
Juan, Argentina, and visits in August 2011 to the Laboratory of High Energy Astrophysics, 
NASA Goddard Space Flight Center, Greenbelt, Maryland, USA, and in February-March 
2012 to the European Space Astronomy Centre (ESAC), Madrid, Spain, the latter funded 
by the COSPAR Fellowship Program. Warm hospitality and tutelage at these institutions 
are kindly acknowledged, in particular from Andy Pollock and Carlos Gabriel at ESAC. We 
would also like to thank Michael Nowak from MIT for useful discussions relating the pileup 
model and the ISIS implementation. 
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Fig. 1. — Light curves of the four observations, ObsID 6615, 6616, 6617, and 6618, in units 
of counts/s. The approximately constant average counts/s that prevails is an indication of 
a low variability degree. 
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Fig. 2. — Chandra MEG spectra of the X-ray binary XTE J1817-330 simultaneously fitted 
in the 11-25 A region using the TBnew model, (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 
6617. (d) ObsID 6618. See Table 2 for further details pertaining these observations. In each 
panel, the black data points are the observations while the solid red lines correspond to the 
best-fit models. Residuals are shown as histograms for each case in units of y 2 . Note the 
large residuals in the oxygen K-edge region, particularly in the neighborhood of the O I and 
O II K a lines at ~ 23.50 A and ~ 23.35 A, respectively. 
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Fig. 3. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a simple functional model (power law and several Gaussian 
profiles), (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. The observed 
absorption features are interpreted as: the neutral O K-edge (~ 23.1 A); the Ka, K/3, and 
Ky lines of O I (~ 23.50 A, « 22.88 A, and ~ 22.68 A, respectively); the Kct, K /3, and Ky 
lines of O II (~ 23.35 A, « 22.28 A, and ~ 22.10 A, respectively); and the K a lines of O III, 
O VI, and O VII (~ 23.10 A, « 22.02 A, and ~ 21.58 A, respectively). See Table 4 for more 


details. 
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Fig. 4. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a powerlaw*warmabs physical model, (a) ObsID 6615. 
(b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 
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Fig. 5. — Comparison of the theoretical photoabsorption cross sections for O 1 (red), O If 
(green), and O Iff (blue) computed by Garcia et ah (2005) which are implemented in the 
warraabs model. The black solid line is the laboratory measurement by Stolte et al. (1997) 
showing peaks of O 1 and O 2 . The vertical dashed lines are placed at the wavelengths of the 
observed absorption lines (see Table 4). Panel (a) displays the original cross sections while 
panel (b) shows the same curves after the wavelength shifts are applied to both O 1 and O II 
(see text for details). 
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Fig. 6. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a powerlaw*warmabs physical model, after the wavelength 
corrections in the O I and O II atomic cross sections are introduced, (a) ObsID 6615. (b) 
ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 


Table 1. Comparison of EBIT and theoretical wavelengths (A) for K lines 


Ion Lower level (J) Upper level (J') EBIT Theory 


O VI 

1s 2 2s(1/2) 

ls2s2p(l/2, 3/2) 

22.0194 ±0.0016“ 
22.374 ±0.008 6 

22.00 c 

22.05 d 

22.03 e 

O V 

1s 2 2s 2 (0) 

ls2s 2 2p(l) 

22.374 ± 0.003“ 
22.370 ±0.010 b 

22.33“ 

22.35 d 

22.37“ 

O IV 

ls 2 2s 2 2p(l/2, 3/2) 

ls2s 2 2p 2 (l/2, 3/2) 

22.741 ±0.005 6 

22.78“ 

22.73 d 

22.75“ 

O III 

ls 2 2s 2 2p 2 (l, 2) 

ls2s 2 2p 3 (l) 

22.071 ± 0.006 6 

23.08“ 

23.05 d 

23.07“ 


a EBIT measurement (Schmidt et al. 2004) 
b EBIT measurement (Gu et ah 2005) 
c HULLAC calculation (Behar & Netzer 2002) 
d i?-matrix calculation (Pradhan et ah 2003) 


e HFR calculation (Garcia et ah 2005) 
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Table 2. Chandra observations used in this paper 


ObsID 

Date 

Exposure (ks) 

Instrument 

Grating 

Label 

6615 

2006 Feb 13 

18 

HETG-ACIS 

MEG 

Observation 1 

6616 

2006 Feb 24 

29 

HETG-ACIS 

MEG 

Observation 2 

6617 

2006 Mar 15 

47 

HETG-ACIS 

MEG 

Observation 3 

6618 

2006 May 22 

51 

HETG-ACIS 

MEG 

Observation 4 
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Table 3. Broadband simultaneous fit parameters 


Model 


Parameter 


Value 


TBNew 

N n (10 21 cm" 2 ) 

1 66"*” 003 
i.OO_ 0 04 

TBNew 

Neon abundance (ANe) 

2.258 

± 

0.080 

TBNew 

Iron abundance (Ap e ) 

1.302 

± 

0.050 

TBNew 

Oxygen abundance (Ao) 

1.178 

± 

0.022 

Power-law 

Normalization 

9.416 

± 

0.064' 



7.039 

± 

0.041 1 



4.613 

± 

0.042' 



1.481 

± 

0.001' 

Power- law 

Photon Index 

1.545 

± 

0.030' 



1.618 

± 

0.030 1 



1.735 

± 

0.029' 



2.242 

± 

0.038' 

SimpleGpile2 

P 

0.050 

± 

0.001' 



0.051 

± 

0.001 1 



0.050 

± 

0.001' 



0.049 

± 

0.001' 

Reduced chi-square 

x 2 

1.208 




Note. — Abundances relative to the solar values of Wilms et al. 

( 2000 ). 

a ObsID 6615 

b ObsID 6616 

c ObsID 6617 

d ObsID 6618 
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Table 4. Absorption line assignments 


Ion 


Transition 

A a (A) 

A b (A) 

0 VII 

Is 2 

4 S - [ls]2p 4 P° 

21.589 ±0.003 

21.60 ±0.06 

0 VI 

2s 2 S - [ls]2s2p 2 P° 

22.022 ±0.003 


0 II 

2p 3 

4 S° - [ls]2p 3 4p 4 P 

22.101 ±0.005 


Oil 

2p 3 

4 S° - [ls]2p 3 3p 4 P 

22.280 ±0.003 


01 

2p 4 

3 P - [ls]2p 4 ( 2 P)4p 3 P° 

22.609 ±0.007 


01 

2p 4 

3 P - [ls]2p 4 ( 2 P)3p 3 P° 

22.686 ± 0.004 


01 

2p 4 

3 P - [ls]2p 4 ( 4 P)4p 3 P° 

22.790 ±0.001 


01 

2p 4 

3 P - [ls]2p 4 ( 4 P)3p 3 P° 

22.884 ±0.004 

22.91 ±0.03 

0 III 

2p 2 

3 P - [ls]2p 3 3 P° 

23.054 ±0.001 

23.13 ±0.09 

0 III 

2p 2 

3 P - [ls]2p 3 3 D°, 3 S° 

23.104 ±0.005 

23.13 ±0.09 

Oil 

2p 3 

4 S° - [ls]2p 4 4 P 

23.358 ±0.002 

23.35 ±0.03 

01 

2p 4 

3 P - [ls]2p 5 3 P° 

23.502 ±0.001 

23.52 ±0.02 


a Present work 


b Sala et al. (2007) 
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Table 5. warmabs simultaneous fit parameters 


Model 

Parameter 

Value 

Warmabs 

JV H (10 21 cm- 2 ) 

1.38 ±0.01 

Warmabs 

Log ionization parameter (log £) 

-2.699 ±0.023 

Warmabs 

Oxygen abundance (Ao) 

o.689t°;°i® 

Power-law 

Normalization 

2.384“ 



2.544 6 



1.511“ 



0.685 d 

Power-law 

Photon Index 

3.044“ 



2.479 b 



2.720“ 



2.617 d 

Reduced chi-square 

x 2 

1.245 


Note. — Oxygen abundance 
Grevesse & Sauval (1998). 

“ObsID 6615 

b ObsID 6616 

c ObsID 6617 

d ObsID 6618 


relative to the solar value of 
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Table 6. Hydrogen column density comparison 


Method 

N h (10 21 cm- 2 ) 

TBnew fit“ 

1 66~*”°' 03 
i.oo_ 0 04 

warmabs fit a 

1.38 ±0.1 

21 cm survey b 

1.58 

21 cm survey 0 

1.39 

Spectral model fit d 

1.55 ±0.05 

Spectral model fit 6 

1.2 

Spectral model fit^ 

0.88 - 0.97 

Spectral model fit^ 

1.0 - 3.0 

Spectral model fit^ 

1.08 ±0.05 


a Present work 
b Dickey & Lockman (1990) 
c Kalberla et al. (2005) 
d Sala et al. (2007) 

e Rykoff et al. (2007), Gierlinski et al. 
(2008), Roy et al. (2011) 

f Millcr et al. (2006a) 

g Torres et al. (2006) 

h Gierlinski et al. (2009) 


40 


Table 7. Equivalent widths comparison for the oxygen Ka transitions in the ISM 


Source 

O I Ka 

(mA) 

O II Ka 

(mA) 

O III Ka 

(mA) 

O VI Ka 

(mA) 

O VII Ka 

(mA) 

XTE J1817-330 “ 

51 + 5 

52 + 5 

16 + 5 

14 + 5 

54 + 4 

XTE J1817-330 6 

70 + 20 

40 + 20 

30 + 30 


40 + 30 

Cyg X-2 c 

67+11 

26+17 

9+8 



Cyg X-2 d 

50 + 3 

28 + 6 


11 + 2 


4U 1543-624 c 

68+11 

41+11 

18+5 



4U 1820-30 c 

70+20 

41+17 

20+10 



4U 1735-44 c 

80+30 

40+20 

24+10 



GX 9+9 c 

80+30 

50+19 

27+13 



4U 1636-53 c 

90+30 

60+20 

33+13 



Cyg X-l (ObsID 3407) c 

95+17 

20+11 

35+15 




a Present work 
b Sala et al. (2007) 
c Juett et al. (2004) 
d Yao et al. (2009) 
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Tabic 8. Oxygen column density comparison 


Source 

O I 

(10 17 cm- 2 ) 

O II 

(10 17 cm- 2 ) 

O III 

(10 17 cm" 2 ) 

O VI 

(10 15 cm“ 2 ) 

O VII 

(10 15 cm" 2 ) 

XTE J1817-330 “ 

3.18 ±0.79 

2.18 ±0.54 

0.35 ± 0.08 

5.06 ± 2.53 d 

45 ± 18 d 





5.99 ± 1.9 e 

988 ± 405 e 





71.12 ± 29.15 f 

4484 ± 1838 f 

XTE J1817-330 6 

5.85 ± 1.75 

0.49 ±0.12 

0.07 ±0.02 



PG 1302-102 c 




0.16 ±0.03 


Mrk 1383 c 




0.38 ±0.07 


ESO 141-G55 c 




0.31 ±0.06 


PKS 2005-489 c 




0.60 ±0.01 



a Calculated from the EWs and the curves of growth (see section 7 ). 

b Derivated from the warmabs model fit using the solar value of Grevesse &; Sauval (1998) (see section 

5 ) 

c Obtained from FUSE data by Savage et al. (2003) 
d Using a velocity dispersion v = 200 km/s 
e Using a velocity dispersion v = 100 km/s 


f Using a velocity dispersion v = 20 km/i 


